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■ Abstract 

o : 

■ Within the nonrelativistic quantum chromodynamics (NRQCD) factorization formalism, we com- 
. pute the relativistic corrections to the exclusive decays of bottomonia with even charge conjugation 

o ; 

T-H , parity into S-wave charmonium pairs at leading order in the strong coupling constant. Relativistic 

^ ' corrections are resummed for a class of color-singlet contributions to all orders in the charm-quark 

^ ■ velocity Vc in the charmonium rest frame. Almost every process that we consider in this work has 

H : 

. 5t , negative relativistic corrections ranging from —20 to —35 %. Among the various processes, the rel- 

ativistic corrections of the next-to- leading order in Vc to the decay rate for Xb2 ric{mS) -\-r]c{nS) 
with m, n = 1 or 2 are very large. In every case, the resummation of the relativistic corrections en- 
hances the rate in comparison with the next-to-leading-order results. We compare our results with 
available predictions based on the NRQCD factorization formalism. The NRQCD predictions are 
significantly smaller than those based on the light-cone formalism by 1 or 2 orders of magnitude. 

PACS numbers: 12.38.-t, 12.38. Bx, 13.20.Gd 
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I. INTRODUCTION 



Among various bottomonium states if, the rjb and Xbj mesons with J^^ = and l"*""*", 
respectively, have a common feature that the charge conjugation parity C is even. Because 
the quantum chromodynamics (QCD) preserves C, these mesons can decay into a pair of 
charmonia hi + h2 with the same C parity. Therefore, possible decay modes of the mesons 
are 

Vb or Xbj -> ip{mS) + ipinS), (la) 
Vb or Xbj VcimS) + rj^inS), (lb) 
Vb or Xbj Xcj + XcJ', (Ic) 

where the spin-triplet S'-wave state ip{nS) is a J^*" = 1 eigenstate. However, some of these 
decay modes like i]}, or Xbi VcijnS) + rjdnS) are forbidden due to the parity conservation 
in QCD. 

Previous theoretical studies include the calculation of the decay rate for the process 
rji, ^ J + J /ijj [1] within the nonrelativistic QCD (NRQCD) factorization formalism j2|. 
This process was proposed as a candidate for the discovery mode of the meson . This 
work was followed by a study of the final-state interaction in Ref. and the next-to-leading- 
order (NLO) QCD corrections in Refs. j^,^]. In the case of the spin-triplet P-wave decay, the 
process XbJ J/ip + J/ip was investigated in Ref. [3] within the light-cone (LC) formalism. 
The authors have extended their LC predictions to various channels in Ref. {s], where they 
also provided the NRQCD predictions at leading order (LO) in the charm-quark velocity Vc 
in the charmonium rest frame. According to the results in Ref. [8|], the LC predictions are 
greater than the NRQCD counterparts although they are in agreement within errors that 
are significant. One of the motivations of this work is to investigate if these discrepancies 
are reduced under relativistic corrections.^ 

In the NRQCD factorization formula, the decay rate is expanded in powers of the velocity 
vq of the heavy quark Q in the quarkonium rest frame. In the case of the LC formalism, the 
amplitude is expanded in powers of the inverse of the hard scale. In the limit of rrib ^ rric 



^ Very recently, the authors of Ref. [9| have reported an NRQCD prediction for the decay XbJ J/ip + J/ip 
including the relativistic corrections of relative order w^, where they used a velocity-expansion scheme that 
is different from the standard NRQCD approach employed in this work. 
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as well as 1 ^ f ^ ^ f^, we can guess that the leading-twist LC prediction and the NRQCD 
prediction at LO in vq are roughly consistent with each other. However, if there are large 
relativistic or QCD corrections, such a naive estimate may fail. For example, the cross 
section for the exclusive production process e^e'^ — )■ J/tp + rjc the B factories suffers large 
relativistic [l^ and QCD [11] corrections. The NRQCD prediction including relativistic 



and QCD corrections is shown to be consistent with the empirical value within errors 10 |. 



One may guess that LC predictions as in Ref. [12|] provide a reasonable answer without 
further corrections. However, as is shown in Ref. [13], the LC calculation contains short- 
distance contributions, which are treated in the context of model LC distributions. These 
short-distance contributions appear in corrections of order (and higher) in the NRQCD 
approach and can be computed from first principles in that approach. Therefore, it is very 
interesting to see what happens in the case of the bottomonium decay into charmonium 
pairs that have similar features as the exclusive process e+e^ J/ip + rjc- 

In this work, we compute the decay rates for the C-even bottomonia rjb and XbJ into 
pairs of S-wave charmonia within the NRQCD factorization formalism. The computation 
is carried out within the color-singlet mechanism of NRQCD because, in these exclusive 
processes, the color-octet channel enters only if at least two hadrons involve color-octet 
contributions that are suppressed.^ In addition, we also neglect the electromagnetic decay 
mode such as r/f, — ?■ 7*7* — ?■ J /ip + J /ip that is tiny compared with the QCD mode.^ 
Relativistic corrections are computed by making use of the generalized Gremm-Kapustin 



relation [18|| that has been employed to resolve the e"'"e^ — )■ J/ip + rjc puzzle [10|. The 
method enables us to resum a class of color-singlet contributions to all orders in Vq. Our 
calculation reveals that the discrepancies between the NRQCD prediction and the prediction 
based on the LC formalism become more severe, particularly in the Xb2 decay into an S-wave 
spin-singlet charmonium pair. This is in contrast to the case of e~^e~ J/ip + rjc. 

This paper is organized as follows. We first describe the NRQCD factorization formula 



^ See Ref. [8| for further discussion regardiiigJJie suppression of the color-octet contributions. 

^ In the case of e+e" — >■ J/ip + J ftp [iJ, uM, the LO cross section is comparable to that Ig, ll7| of 
e^e^ J/ijj + rjc because of the enhancements due to the photon fragmentation and large collinear 
emissions from the electron line in the forward region. In the case of rjb — >■ 7*7* — > J ftp + J ftp only the 
photon fragmentation contributes and the collinear enhancement is missing because of the large bottom- 
quark mass. In addition, the fractional electric charge of the bottom quark makes the rate insignificant. 
See also Ref Q. 
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for the exclusive C-even bottomonium decay into an S'-wave charmonium pair in Sec. |TT1 In 
Sec. [TTll we present the strategy to compute the short- distance coefficients for the NRQCD 
factorization formula. The analytic results for the decay rates at LO and the corrections of 
NLO in v1 are given in Sec. IIVI Our final numerical results for the decay rates, in which a 
class of the color-singlet contributions are resummed to all orders in f ^, are listed in Sec. |V] 
and compared with available predictions. Finally, we summarize the work in Sec. IVII 



II. NRQCD FACTORIZATION FORMULA 

The exclusive decay of a bottomonium into a pair of charmonium states involves the 
annihilation of a hb pair followed by the creation of two pairs of cc. One may guess that the 
generalization of the NRQCD factorization [2] for the electromagnetic decay or light-hadronic 
decay into this exclusive mode is possible.^ If we assume that the NRQCD factorization is 
valid for the exclusive decay of a bottomonium H into a charmonium pair hi + h2, then 
the decay rate can be expressed as a linear combination of the products of nonperturbative 
NRQCD matrix elements with numerous spectroscopic states. According to the velocity- 
scaling rules of NRQCD [2I, these matrix elements are classified in powers of vq. These 
exclusive decay modes are dominated by the color-singlet channels as is stated in the previous 
section and we restrict ourselves to the color-singlet contributions. The typical velocity of 
the bottom quark Vb in the initial state is significantly smaller than those of the final-state 
charmonia, ~ 0.1 ^ f ^ ~ 0.3 so that we neglect the relativistic effects of the bottom 
quark while we include the relativistic corrections of the charm quarks in the final-state 
charmonia. 

Within the color-singlet mechanism at LO in vq, there are NRQCD matrix elements 
{H\Oi\H) for the bottomonium decay and (0|(9j^'|0) for the charmonium production that 
involve the decay rate r[H ^ hi + for z = 1 and 2. The spectroscopic states for the four- 
quark operators Oi and are identical to those of the corresponding hadrons.^ In order 
to describe the relativistic corrections to the charmonium state, we denote (OlC^'*^^, „.)|0) 
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For example, NRQCD factorization theorems for the exclusive quarkonium productions in e+e annihi- 



lation and B decay have been proved 19l42l|. 



^ For the initial-state bottomonia '^^'^^Lj = ^Sq and ^Pj with J = 0, 1, or 2 for H = r]b and Xbj, 
respectively. In the charmonium case, ^^^^Lj — ^So and ^Si for hi — r]c and J/ip, respectively. 
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as the NRQCD matrix element for the charmonium production that is of relative order 
in comparison with the LO matrix element (O|0^'|O). The spectroscopic state of 



V, 



2{mi+ni) 



the four-quark operator 0^^^^_ is again identical to that of hi. As a result, the NRQCD 
factorization formula for V[H ^ hi + hQ\ can be expressed as 

V[H + h,] = {H\Oi\H) J2 C(™i,nO,(n^.,n.)(0|Ot;„^,„^)|0)(0|0?^^(^^^„^)|0), (2) 

mi,ni 

where C(^mi,ni),{m2,n2) is the short-distance coefficient which is insensitive to the long-distance 
nature of the hadrons H, hi, and h2. These factors can be computed perturbatively in 
powers of the strong coupling a^. 

The LO color-singlet NRQCD four-quark operator Oi for the annihilation decay of the 
heavy quarkonium with the spectroscopic state "^^^^Lj is of the form 

Oi = ^^}CC'^'Lj)xx^lCC'^'Lj)i^, (3) 

where ip and are the Pauli spinor fields that annihilate Q and Q, respectively. For 
a bottomonium (charmonium), it is understood to he Q = b (c). Here, the operators 
^^25+1^ j) are defined by 

ICCSo) = 1, (4a) 
IC'CSi) = a\ (4b) 
/C(3Po) = ^(-t'^-^), (4c) 

rCPi) = ^(-t"^ X ay, (4d) 

V^CP2) = -i^(V^), (4e) 

where 1 is the identity matrix for the spin and color space, cr* is the Pauli matrix and D 
is the gauge-covariant derivative. The notation A^^^^ in Eq. f He|) represents the symmetric 
traceless component |(A*-' + A^^) — of a Cartesian tensor A*-^ . The NRQCD matrix 

element {H\Oi\H) for the decay is averaged over the spin states of H. 

The NRQCD four-quark operators C^i(m„) in Eq. ([2]) for the S-wave charmonium pro- 
duction can be expressed as 



O'um^n) = \ X^/C™('^+^L,)^(5^aI^a,,)^X(2^+^L,)x + H.c. 



(5) 
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where /C„ are defined by 

ICnCSo) = (-fV)^", (6a) 
ICnCSi) = (6b) 



As is explained earlier, the Pauli spinor fields are for the charm quark. We have listed 
only a class of the operators that contain ordinary derivatives rather than covariant deriva- 
tives. The neglect of the operators with gauge fields contributes first at relative order 
in the Coulomb gauge in which the matrix elements are evaluated.^ Applying the vacuum- 
saturation approximation, we can simplify these NRQCD matrix elements as 

(O|0te^,„)io) = {2j+i){Q\x^iCmmmx^iCnm* + o{vt), (?) 

where J is the total angular momentum of the charmonium h. The spin-multiplicity factor 
2 J + 1 appears because the spin states of the produced hadron h are summed over in Eq. 

In carrying out the resummation of the relativistic corrections, it is convenient to define 
the ratio {q^"')h of the NRQCD matrix element of relative order f^" to the LO matrix element 

as 



This ratio is independent of the polarization of h. In Ref. 



22| a generalized version of the 



Gremm-Kapustin relation 2J] was derived: 

(9^"). = {q'K- (9) 

This relation holds for the matrix elements in spin-independent-potential models. Thus this 
relation holds for both spin-singlet and -triplet states independently of the index i up to 
corrections of Vq that break the heavy-quark spin symmetry. This relation has been applied 
to determine the NRQCD matrix elements for the S'-wave quarkonium states precisely 



25| and to resum the relativistic corrections in various quarkonium processes [10|, |23|, |26 |. 
Although there is an intrinsic limitation that the resummation of relativistic corrections with 
only the QQ Fock-state contributions eventually has the predictive power up to corrections 
of relative order Vq [isj, the method is still useful to improve the convergence in a process 
involving significant relativistic corrections. Such an example is the exclusive J/ip + tjc 



6 See Refs. 



22 



23| for further discussion. 
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production in e~^e~ annihilation 10|. We shall find in Sec. |V] that a large correction is 
indeed observed in the process Xb2 Vc + Vc- Because of the large relativistic corrections of 
relative order v"^, the theoretical prediction for the cross section can even be negative. This 
problem will be cured by resumming the relativistic corrections. 

By employing the vacuum-saturation approximation and the generalized Gremm- 
Kapustin relation (Q, we can express the NRQCD matrix element for the charmonium 
production as 

{0\Ol(rn,n)\0) = i2J+l){qX^-{h\O^\h)+O{vt), (10) 

where the errors of order are from the vacuum-saturation approximation. Then the 
NRQCD factorization formula ([2]) for the decay is simplified into the form 

T[H + h2] = {H\0,\H){h\Oi\h){h2\Oi\h2) dn,,nAQX'MXl^ (11) 

ni,n2 

where we have redefined the short- distance coefficient as (i„i,n2- The state \H) in the NRQCD 
matrix elements has the nonrelativistic normalization {H{P)\H(P')) = {2tt)^ 6^^\P — P'). 

The short- distance coefficients dni,n2 ^i'g determined by the perturbative matching. In 
fact, we can construct the NRQCD factorization formula for the amphtude AH^hi+h2 under 
the vacuum-saturation approximation. Once we replace the hadrons H, hi, and /i2 with the 
perturbative heavy-quark- antiquark states bb, cci, and CC2 with the same spectroscopic states 
as those of the corresponding hadrons, respectively, then we find that the QQ counterpart 
Ai£^cci+cc2 is calculable perturbatively as 



X 

ni,ra2 



A6^,5,+ec-2 = V(^^l^ll^^)(cCl|(^l|cCl)(cC2|Ol|cC2) J] a„,,„,g^lqf % (12b) 



ni,n2 



where a„j,„2 is the short- distance coefficient at the amplitude level and qi is half the relative 
momentum of the i-th cc pair. While the hadron states like and \hij are normal- 
ized nonrelativistically, we use the relativistic normalization for the heavy-quark state \Q)'- 
{Q{p)\Q{p')) = (2vr)^2(mQ + p'^Y^'^5'^^\p — p'). The normalization factor y/8mHmh^mh2 
was introduced to make the amplitude AH^hi+h2 have the relativistic normalization like 
Af,i^cci+cc2- Here, mn and m/j. are the masses of H and hi, respectively. 
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Because the short-distance coefficients are insensitive to the long-distance nature of the 
hadrons, the factor ani,n2 in Eq. f lT2|) must be common to both Af,i^cci+cc2 AH^hi+h2- 
The QQ NRQCD matrix elements in Eq. ( ]12b|) are calculable perturbatively: 



|Ci|66) = 2iVe(2E)2, (13a) 
{cci\0i\cci) =2N,{2Eif, (13b) 

where E and Ei are the energies of the quark in the QQ rest frame for the bb and cq 
pairs, respectively. In such a way, one can determine the short- distance coefficients am.ng 
and dni,n2- the normalization of the P-wave matrix element (fT3|) we have suppressed the 
factor |qp that comes from the derivative operator.^ 



III. PERTURBATIVE MATCHING 



In this section, we present the method to compute the short- distance coefficients (i„i,„2 
the NRQCD factorization formula for the decay — ?■ /ii + /12 by employing the perturbative 
matching onto the full-QCD amplitude for bb — )■ cci + CC2- 



A. Kinematics 



The momenta for the quarkonia H, hi, and h2 are chosen to be P, Pi, and P2, respectively. 
We list the definitions of the variables for the bottomonium H or the bb pair without any 
index. Variables with the index i = 1 or 2 indicate that they correspond to the charmonium 
hi or the ccj pair. For convenience we list the definitions for various variables for a QQ pair 
system without any index unless it is necessary. 

We denote P and q by the total and half the relative momentum of a QQ pair, respectively. 
Then the momenta for the quark (p) and the antiquark (p) are expressed as 

p=^P + q, (14a) 
p=^P-q. (14b) 

^ See, for example, Eq. (60) of Ref. 
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It is obvious, in the rest frame of the QQ pair, that P and q are orthogonah P ■ g = 0. In 
that frame. 



P = (2E,0), (15a) 

g = (0,q), (15b) 

p={E,q), (15c) 

P={E,-q), (15d) 

where E = ■\JfnQ + is the energy of the Q or Q in the QQ rest frame and mg is the mass 
of the heavy quark Q. We assume that the Q and Q are on their mass shells so that 

p'^ = = ml, (16a) 

p2 = 4El (16b) 



In addition, we define the ratio 



which is useful to simplify expressions. 



r = — , (17) 

ml 



B. Spin and color projectors 

The Feynman diagrams for the process bb — t- cci + CC2 at LO in are shown in Fig. [TJ 
The corresponding perturbative amplitude for the process bb — )■ cci + CC2 is of the form 




FIG. 1: The Feynman diagrams for the exclusive decay of a C-even bottomonium decay into a pair 
of charmonium states at LO in as- 
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■ 4 

A = , _ .3/ , _ ,^ v{p, s)13f,^u{p, s)u{pi, Si)C''v{p2, S2)u{p2, S2)C''v{pi, Si), (18) 

[Pl +P2)^{P2+Pl)^ 

where Qs is the strong couphng, the denominator factor is from the gluon propagators and 
s and s (sj and Sj) are the spins of b and b (c and c in the i-th cc pair), respectively. B'^'' 
contains the bottom-quark propagator and the gluon vertices to the bottom-quark line. 
is the gluon vertex to the charm-quark pair cc and we have suppressed the color indices. 
Both B^'^ and act on spinors with both Dirac and color indices. 

In order to compute the amplitude for the perturbative process bb — )■ cci + CC2 with the 
appropriate spectroscopic states, it is convenient to use the projection operators. In this 
work, we consider only the color-singlet contributions that can be projected out by replacing 
the color component of the outer product of the spinors for Q and Q in each QQ pair with 
the color-singlet projector 

where 1 is the unit color matrix. 

The spin-singlet and -triplet components of each QQ state can be projected out by making 
use of the spin projectors. After multiplying corresponding Clebsch-Gordan coefficients to 
the spin component of the outer product of the spinors for each QQ pair, one can find the 
spin-singlet and -triplet projectors Hi and II3 for the QQ decay and Hi and TI3 for the 
QQ production, respectively. The spin projectors that are vahd to all orders in the relative 



momentum can be found in Refs. 



20 



28|. 



Hi = -iV(|^ + mQ)(f +2E)75(^-mQ), (20a) 

n3a(A)e-(A)= N{]l, + mQ)ir+2E)^i\){p-mQ), (20b) 

ni= N{p-mQ)^,i^+2E){^ + mQ), (20c) 

n3„(A)e*°(A) = Nip-mQ)fi\)i^+2E)iP + mQ), (20d) 

where the normalization factor is 

N = ^ , (21) 

AV2E{E + mQ) 

if we choose the relativistic normalization for the spinors. e(A) is the polarization four-vector 
of the spin-triplet state with the helicity A. 
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C. Projection of S- and P-wave contributions 



We can extract the color-singlet amplitude Ai with appropriate spin states of the three 
quarkonium states from the amplitude A as 

• 4 

^1 = ^— -w^— -^Tr[s^.n|,5] TTic-uurnu], (22) 
[pi +P2r{p2+pir 

where and are those in Eq. ( |T8|) and 11 (11) is the direct product of the color and 
spin projectors defined in Eqs. f|T9|) and f l20l) : 

n = TTi ® (Hi or Hg^e") , (23a) 
n = TTi ® (Hi or Ilaae*") . (23b) 

The trace in Eq. f l22p is over the color and spin indices. 

As the next step, we need to pull out the L-wave amplitude, where L = S or P, from 
the color-singlet amplitude Ai in Eq. f l22|) with correct spin states. In the case of the initial 
states H = rjh and Xbj, we need to project out the spin-singlet S'-wave and spin-triplet P- 
wave states, respectively. As we have stated earlier, we consider only the contributions of 
LO in Vb in the bottomonium sector. Because the dependence of the momenta P and q for 
the bb pair is isolated in the tensor B'^'^ in Eq. f l22l) . the projection of the bb state can be 
made only with this factor as 



(24a) 

q=0 



C.(3P.)=^f|^Tr[i3-(vr,®n3,|,,)] 



(24b) 

q=0 



where Vj'^ projects the total angular momentum state J in the P-wave spin-triplet contri- 
bution that are defined by 

Vf = i=r^ (25a) 

Vf = 1^^''""' PM>^1 (25b) 
= e"^(A). (25c) 

Here, e'^(A) and e"^(A) are the polarization vector and tensor for the ^Pi and states with 
the hehcity A, respectively. In the normalization of B^^^^j,p^^ in Eq. fl24l) . we have suppressed 
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the factor |qr| to make it to be consistent with the normahzation of the matrix element (fT3l) . 
The tensor is defined by 

pa pl3 

= ^g»^ + __. (26) 

Finally, we summarize the method to project out the S-wave component of the amplitude 
from Ai in Eq. f l22p for the two-charmonium final states that are either a spin-singlet or a 
triplet including relativistic corrections. Because Ai in Eq. (!22l) is the color-singlet amplitude 
of the ccj pairs with appropriate spin states, we only need to pull out the S-wave component 
that is independent of the direction of Qi but may depend on qf. We notice that the trace 
factor that includes depends on both qi and q2 through the bottom-quark propagator. 
The trace factor for the charm-quark pairs and the factor of the gluon propagator l/[{pi + 
V2Y{p2 + PiY] have the dependence on both qi and q2. After taking the average over the 
directions of qi and ^2, we find our final expression for the perturbative amplitude Ahi^cci+cc2 
as 



' dVLidVL2 Tr[C^n|cg2Cn| 
(47r)2 ^"^i(pi+p2)2(p2+Pi) 
where B'^^^ is defined in Eq. fl2^ . Note the hb and cq pairs in Eq. fl27|) have definite 

spectroscopic states that are suppressed. Here, dVti is the solid-angle element that represents 

the direction of qi. 

D. Short-distance coefficients 

The perturbative amplitude ^bb-s>cci+cc2 depends only on m^,, mc, q]^ and polarizations A 
and Aj of hh and ccj pairs, respectively. By making use of the normalization of the NRQCD 
matrix elements for the QQ states in Eq. ( !T3|) . we can find the NRQCD factorization formula 
for the amplitude of the decay H{\) — )■ /ii(Ai) + h2{X2) 



{H\Oi\H){h\0,\h){h2\Oi\h2) 



where we have suppressed the helicities of the hadrons that are the same as those for the 
QQ pairs. Note that qf in the perturbative amplitude Ai)i^cci+cc2 is replaced with the ratio 

Now we finally find the NRQCD factorization formula for r[H ^ hi + h2]: 

T[H ^hi + h2] = j ^ ^^^mg 5Z \^H{X)^h^{\^)+h2{\2)? , (29) 
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where J is the total angular momentum of H , J d^2 is the phase space of the final state and 
the summation is over the polarizations of the initial and final states. The factor 2m// in the 
denominator cancels that of the squared amplitude.^ After summing over the polarization 
states of the particles, the squared amplitude becomes invariant under rotation and the 
phase-space integral becomes trivial 

where S is the symmetry factor, S* = 1 (S* = 2) for hi ^ {h\ = h2). Pcm is the magnitude 
of the three-momentum of hi in the H rest frame 



where A(x, y, z) = + + — 2{xy + yz + xz). Note that we use the physical masses m// 
and m/i- in evaluating the phase space $2- This prescription respects the physical endpoints 
of the phase space without spoiling the gauge invariance. 

In evaluating the S'-wave amplitude f l27|) , we follow the strategy given in Ref . [l^ . In order 
to investigate the convergence of the series in f^", we also provide the fixed-order prediction 
at LO and NLO in v^. To compute the fixed-order relativistic corrections, we first expand 
the amplitude ( l27ll in powers of Qi and then take the angle average of gj-dependent tensors 
by making use of the following formulas: 



\^ = 0, (32a) 



47r 



^gf^r = |r , (32b) 

"^^^^^.— O, (32c) 



47r 

..u .V B Q 



,4 



4^ q^q^qTq? = ^ ( irir + irc + irir ) , (32d) 

where /f^ is the same as I^^ defined in Eq. fl26|) except that we replace P with Pj. 

Our final expression for the short- distance coefficient (i„j.„2 in the NRQCD factorization 
formula flTT]) for the decay rate V[H ^ hi + h2] is 



dm ,712 



niW d^d^'- 



/A,Ai,A2 \'^bb{\) — >cci(Ai)+cc2(A2) I 

(2J+l)(2iV,)3(2mb)2(2Ei)2(2E2)^ 



• (33) 



See Eq. ^2M 
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In the summation of the polarization states for the spin-1 and -2 quarkonia, we use the 
following formulas: 

^e°(A)e*^(A) = J"'^, (34a) 

A 

J2 e'^^ixy^'^W = -(r^/^^ + r'^i^'') - (34b) 

A ^ ^ 

where I'^" is defined in Eq. ( I26l) . For the spin-1 and -2 charmonium hi, I^" must be replaced 
with If". 

IV. THE DECAY RATE UP TO NLO IN vl 

Now it is straightforward to compute the decay rate based on the strategy and techniques 
described in Sec. IIIII In this section, we present the analytic expressions for the decay rates 
V[H ^ hi + /i2] for H = rjh, Xbo, Xb2 and for hi = ip{mS), rjcijnS) up to NLO in v^. In the 
case of Xbi — ^ J /''P + J li^^ the LO contribution first appears at order t>^ and the decay rate is 
of order v^. This process is so strongly suppressed that we present only the LO contribution. 

We write the decay rate of the form 

V[H /ii + /12] = T^'^iH ^hi + /i2](l + R), (35) 

where F^^ is the LO contribution and R is the ratio of the relativistic corrections to the LO 
contribution. We also define R2 as the corresponding ratio at NLO in v"^. In the following, 
we use the symbols i/'i and r^j for hi = 'ip{mS) and ric{mS), respectively. 



A. % Vl + ■02 



Our analytic results for F^'-' and R2 for rjh ^ ipi + ip2 are given by 



40967r3af (1 - 4r)3/2 



65Qlml'^rS 

-1 



15m2r(l-4r)((q%, + (g%J 



{r],\OM{HOi\^i){HOi\H 
^ (36a) 
3((q^)J, + (g%)(3 + r + 8r2) 



(36b) 
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where the symmetry factor S = 2 for ifji = ijj2 and, otherwise, S = 1. The variable r is 
defined in Eq. f|T7|) . For simphcity, we have put nih^ = 2Ei in Eq. fl36|) including the phase- 
space factor $2- However, when we present the numerical results, we will use the expression 
f PT]) to evaluate $2- As shown in Eq. (15^ . r^*-*[?7b ipi + tp2] is of order v^. In the color- 
singlet spin-singlet bb contribution in Eq. (I24ap . the Dirac trace of Tr[i3^,^n|^5] in Eq. ( 12^ 
must be proportional to e^^"^kiak2i3, where ki = pi + p2 and k2 = P2+V1 are the momenta 
for the virtual gluons. In the limit Vc = 0, however, the two momenta become identical to 
each other ki = k2 = {Pi + P2) /2 = {nib, 0, 0, 0), and therefore the amplitude vanishes before 
the vector indices are contracted to the polarization four- vectors of ipi and tp2- As a result, 
the leading contribution to the amplitude begins at order v"^ and r^*-* [T]h — )■ V^i + 1P2] is of 
order f^. The expression for T^'^[rih tpi + 1P2] in Eq- dSn]) agrees with Eq. (22) of Ref. 
The result for R2 is new. 



B. xbo m + m 



The results for Xbo Vi ~^ V2 are given by 



or ^ ^ 1 2(l-3r + 3r^ + 8r3)((gX + (q%) 

While the amplitude for rjf, ^ ipi + ip2 begins at order v"^, the LO contribution to the 
amplitude A^f^^r]i+r)2 begins at relative order v^. The analytic expression for r'"^ was also 
given in Ref. [8|, where the authors evaluated the phase-space factor $2 in the limit of 
fT^hi 0. Aside from this difference in the phase-space factor, the analytic expression for 
T^'^lrji, — 7- r/i + 772] in Eq. fl37j) is smaller than that in Ref. by a factor of 2. 



C. Xb2 ^ V1+V2 

The results for Xb2 Vi + V2 are given by 

r^°[X62->r^i+r72] = io935^io^g {Xb2\0,\xb2){vi\OM{V2\OM, (38a) 

R2[Xb2 ^V.+ V2] = (1 - 4r) " ^''^^ 



15 



As in the case of T^'-'[rih — > r?i + 772], T^^[xb2 Vi + V2] in Eq. ( 138|) is smaller than the 
corresponding result in Ref. [Sj] by a factor of 2. The result for R2 is new. 



D- Xbo ■01 + ■02 



The LO decay rate and relativistic correction factor R2 for Xbo — ^ "^i + ■02 are 



r^°[X6o 01+^2] = 2187ml^rS (x.o|Oi|x.o) 

x(0i|Oi|V'i)(02|Ci|02), (39a) 
^2 Xw ^ ^1 + ^2 = 271 o , OQ 2 — ■ 39b 



The results in Eq. (139|) agree with those Ref. jsl.^ Like T^'^[rii) rji + 772] and r'"*-'[x62 — > 



^1 + ^2], r^^[xw ~> ^01 + ■02] in Eq- ( p9|) is smaller than the corresponding result in Ref. js] 
by a factor of 2. The result for R2 is new. 

E. Xbi ^ ^1 + V'2 

The LO amplitude for Xbi — ^ '0i + 02 is of order to make the decay rate of order v^. 
Because the process is highly suppressed, we only present the decay rate at the LO in Vc- 

r^°[x.i ^01 + ^^2] = 4428675^i8r4g ' (x.i|gilx.i)(V^i|gil^i)(02|0i|02) 

X [9r(8r + 25)i{q%^ + {q%J + 480r(r + 2){q')^^{q')^^ 

X {{q%, + {q')l) + 2(472r2 + 1375r + 1600)(q^)J^(g2)JJ . (40) 

The expression r'"°[xf,i — ^ ^01 + ^2] in Eq. fHOj) is new. 



Since the authors in Ref. ^9'] adopted a different relativistic expansion, the comparison has been made 
after expanding mj/^ = 2-\/m^ + in their expressions in powers of . 
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F. Xb2 -01 + ■02 



The LO decay rate and relativistic correction factor R2 for Xb2 "ipi + 'ip2 are 
T.LOr , ,1 10247r3a^(13 + 56r + 48r2)v/r^4f , , , 

r-[x.2 ^ + H = io935..,^o,^ ^^-"^^"^-^ 

x{^,\Ol\^^J,){^|J2\0,\^P2), (41a) 

^2 X62 ^ ^/'l + ^/'2 = 2 /ofi I o 5^T^ 5^r^\ • (41b) 



While the resuhs in Eq. ( |4TI) are in agreement with those in Ref. j9|, T^'^[xb2 i^i + ^2] is 
smaller than the corresponding result in Ref. jsj by a factor of 2. 



V. RESUMMATION OF RELATIVISTIC CORRECTIONS AND NUMERICAL 
RESULTS 

In Sec. IIVI we have listed the NRQCD factorization formulas for the decay rates r[iJ — >■ 
hi + /i2] for various processes, in which the LO predictions and the NLO corrections with 
respect to are included. In this section, we provide our predictions for these decay 
rates including relativistic corrections to all orders in in which a class of color-singlet 
contributions are resummed. Still we neglect Vb- The approach that was employed in Sec. IIVI 
at fixed orders in cannot be used to resum the relativistic corrections to all orders in f^". 
Instead of carrying out the expansion order by order, we evaluate the average of the 



amplitude 
IV of Ref. 



[27D over the direction of Qi numerically following a previous analysis in Sec. 



10|.^° Because the calculations are carried out for the perturbative amplitude, 
expressions are for the parton process bb — )■ cci + CC2 rather than the process H ^ hi + h2. 

For simplicity, we carry out the calculation in the rest frame of the initial bb pair. In this 
frame the explicit components of the four- momenta for the bb, cci, and CC2 are given by 

P* = (2mb, 0,0,0), (42a) 
Pi* = (El, 0,0, Pcm), (42b) 
P; = (^2, 0,0, -Pcm), (42c) 



10 



In some simple cases, complete analytic expressions that contain the relativistic corrections resummed to 
all orders in are known. See, for example, Refs. 



23 



25 
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where the superscript in a four-vector V* indicates that the four- vector is defined in the bb 
rest frame, Pcu and Ei are the momentum and the energy of the cq pair: 

AV2[(2m,)^(2E0^(2i?2)'] 



-PcM 



Arrib 



We can parametrize the momentum qi in the rest frame of the cq rest frame as 



li = 1 9i I (0 5 sin cos , sin 9i sin 0j , cos 9i 



(43a) 
(43b) 

(44) 



where di and (pi are the polar and azimuthal angles of qi in the cq rest frame. After boosting 
qi and q2 to the bb rest frame, we find that 



Qi = |gi|(+7i/3i cos6'i, sin 6^1 COS01, sin ^1 sin 01, 7i cos6'i), 
?2 = 1 92 1 ( -72/^2 cos 6^2 , sin 62 cos 02 , sin 62 sin 02 , 72 cos 62 ) , 



where 



7i 



-PcM 



(45a) 
(45b) 



(46a) 
(46b) 



In evaluating the angle average of the amplitude, we have to choose the numerical values 
for the input parameters such as the strong coupling as, heavy-quark mass mg, NRQCD 
matrix elements and physical masses mfj and rrih^ of the hadrons. In this work, we take 
«s(^fe) = 0.215 that is, within uncertainties, consistent with previous analyses on Xbj 
CC + X I27I, T(IS') CC + X [29!, and rib — 77 25]. The masses of the heavy quarks c and 



b are chosen to be the one-loop pole masses =1.4 GeV and nib = 4.6 GeV, respectively, 




29|. The most recent fit for the matrix 



following previous NRQCD analyses 
element for rib{lS) can be found in Ref. [25|], which is consistent, within uncertainties, with 
the values used in Ref. Sj. That of the P-wave states Xbj can be found in Ref. 27|. In 
the case of J/ip, the color-singlet NRQCD matrix elements were fit to the electromagnetic 
decay rate of the state in which the relativistic corrections are resummed to all orders in 



,2n 



including the order-a^ corrections. These values can be found in Ref. . For the spin- 
singlet states ric{lS) and ric{2S), we quote the recent values in Ref. 3^ that were fitted to 
the light-hadronic and electromagnetic decay rates with the NRQCD factorization formula 
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TABLE I: The color-singlet NRQCD matrix elements (ME) {h\Oi\h) that is LO in vq and the 
ratio {q^)h for various heavy quarkonia h. We keep only two digits of significant figures. 



ME \ /i 


vb m 


Xbj [27] rj, [30] 


Vc{2S) [30] 


J/^|; [18] 




3.1 GeV^ 


2.0 GeV^ 0.40 GeV^ 


0.20 GeV^ 


0.44 GeV^ 






0.45 GeV^ 


0.50 GeV^ 


0.44 GeV^ 



TABLE II: The decay rate for the process H ^ hi + h2. T is the NRQCD prediction at LO 
in Vh and Vc- F is that in which the relativistic corrections of all orders in f^" are resummed. The 
relativistic correction factor R defined in Eq. (j35p resummed to all orders in v^^ and i?2 at NLO 
order in v^. F^^ represents the LC prediction based on the formula given in Ref.[8]. r^o, r and 
r'"*-' are in units of eV. 



channel 


pLO 


r 


R2 


R 


pLC 


% ^ 


' J/V + J/iJ 


0.58 


0.45 


-27% 


-22% 




Xbo - 




28 


21 


-31% 


-25% 


51 


Xbo - 


-^rjc + Vci2S) 


28 


23 


-33% 


-17% 


76 


Xbo - 


rjc{2S) + r?e(25) 


7.0 


6.4 


-34% 


-10% 


30 


Xb2 - 




0.79 


0.32 


-93% 


-59% 


16 


Xb2 - 


r/c + Vc{2S) 


0.79 


0.35 


-98% 


-56% 


21 


Xb2 - 


Vc{2S) + rjc{2S) 


0.20 


0.09 


-103% 


-54% 


8 


Xbo - 


J/i; + J/il^ 


18 


15 


-20% 


-16% 


79 


Xbi - 


J/i; + J/tlj 


1.0 X 10"^ 


3.1 X 10"'' 




-70% 




Xb2 - 


J/i; + J/tlj 


45 


35 


-34% 


-23% 


270 



with the accuracies of order a^f^- All of these NRQCD matrix elements with corresponding 
references are listed in Table. HI The values for the physical masses for the involving hadrons 
are taken from Ref. Isil as m^, = 9.3909 GeV, m^.^ = 9.85944 GeV, m^^^ = 9.89278 GeV, 
m^,, = 9.91221 GeV, m^^ = 2.9803 GeV, mj/^ = 3.096916 GeV, and m^,(2s) = 3.637 GeV. 

Our NRQCD predictions for the decay rate r[H — t- hi + /12], in which the relativistic 
corrections are resummed to all orders in f^" keeping Vb = 0, are listed in Table. |TT]in units 
of eV. Here, is the prediction at LO in Vb and Vc- In order to demonstrate the significance 
of the resummation, we list the values for R and R2 that are defined in Eq. fl35|) . T^'^ is 
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the prediction based on the LC formula quoted from Ref . [Sj] . The relativistic corrections to 
all of the processes listed in Table. [TTl are negative. According to R2 in Table. Ull the NLO 
relativistic corrections are between —20 and —35 % of the LO prediction except for Xb2 decays 
into spin-singlet S-wave charmonium pairs i]i+ri2 that are more than —90 %. Especially, the 
decay rate for Xb2 '7c(25') +ric{2S) becomes negative at NLO in t>^. The resummed results 
for the decay rates are greater than the NLO-corrected values. The resummation of the 
relativistic corrections to all orders in f^" makes all of the rates T[xb2 Vi + V2] positive. 
The previous predictions of the decay rates T^^ in Ref. jsl based on the LC formula are 
greater than our final results F that include resummation of relativistic corrections to all 
orders in f^".^^ For Xbo Vi ~^ V2, the LC results are greater than ours by factors ranging 
from 2 to 5. Especially, in the case of Xb2 Vi + V2, the LC results are greater than ours 
by factors ranging from 45 to 90. In the case of Xbj — ^ J/ip + J/ip the factors are 5 and 8 
for J = and 2. 



VI. SUMMARY 

We have presented the NRQCD predictions for the decay rates of the bottomonia rib and 
Xbj, that are even eigenstates of the charge conjugation parity C, into S'-wave charmonium 
pairs. The short- distance coefficients for the NRQCD factorization formula are obtained at 
LO in Os- A class of relativistic corrections of the charm quark in the final-state charmonia 
is resummed to all orders in v'^^ by making use of the generalized Gremm-Kapustin relation 



in Refs. 
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22| that is valid in spin-independent potential models and we have neglected 
the motion of the b quark in the initial bottomonium. 

The results show that the relativistic corrections to the decay rates at NLO in are all 
negative. Severe relativistic corrections are observed especially in Xb2 Vi + V2 ioi rji = 77, 
or Tjci^S). The decay rate for Xb2 — ^ '7c(2S') + ?7c(2S') at NLO in vl is even negative. In almost 
every case, the resummation of relativistic corrections of a class of color-singlet contributions 
eventually gives sizable growth of the decay rate in comparison with the NLO predictions. 



c 



Parts of our results are in agreement with previous NRQCD predictions in Refs. 



While we list only the central values of our predictions in Table. |TI1 the authors of Ref. |8| have provided 
the results with uncertainties that are huge. The discrepancies between NRQCD and LC predictions may 
be partially relieved under those uncertainties. 



20 



In comparison with a previous analysis based on the LC formahsm in Ref. js], our NRQCD 
results severely underestimate the decay rates by 1 or 2 orders of magnitude. It is, therefore, 
very important to pin down the source of such significant discrepancies between the NRQCD 
factorization and LC formalisms. We recall previous theoretical studies on the exclusive 
process e^e^ Jjjh + 7]^ aX B factories. The NRQCD predictions of the cross section at 
I, [l^ severely underestimated the empirical values 32-34|. According 



LO in a. and Vc [16| 



to Ref. 



121], the LC prediction, which is much larger than that of NRQCD at LO in as 



and Vci can explain the measured cross section. However, as is shown in Ref. [13], the LC 
calculation contains short-distance contributions, which are treated in the context of model 
LC distributions. These short-distance contributions appear in corrections of order (and 
higher) in the NRQCD approach and can be computed from first principles in that approach. 
We know that the problem has been resolved within the NRQCD factorization formula 



in combination with the NLO corrections in a. 



corrections 



Q withi 



111 and the resummation of relativistic 



within errors. 
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